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Abstract 

We report a quantum magnetotransport signature of a change in Fermi surface topology in the 
Rashba semiconductor BiTel with systematic tuning of the Fermi level Ep. Beyond the quantum 
limit, we observe a marked increase/decrease in electrical resistivity when Ep is above/below 
the Dirac node that we show originates from the Fermi surface topology. This effect represents 
a measurement of the electron distribution on the low-index (n = 0,-1) Landau levels and is 
uniquely enabled by the finite bulk k z dispersion along the c-axis and strong Rashba spin-orbit 
coupling strength of the system. The Dirac node is independently identified by Shubnikov-de Haas 
oscillations as a vanishing Fermi surface cross section at k z = 0. Additionally we find that the 
violation of Kohler’s rule allows a distinct insight into the temperature evolution of the observed 
quantum magnetoresistance effects. 
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Dirac’s matrix equation of relativistic electrons [I] has in recent years been found to be 
profoundly linked to the dynamics of electrons in solids. Application of this formalism has 
proven key in describing the intertwined (pseudo)spin degrees of freedom [2H3]. One striking 
experimental implication of the linear dispersion is that in laboratory magnetic fields, it 
significantly enlarges low index Landau level energy separations compared to parabolic bands 
generated by a comparable tight-binding transfer integral, making emergent Dirac Fermion 
systems an ideal platform to study the quantum Landau level effects. Observations of this 
range from the discovery of Shubnikov-de Haas (SdH) oscillations in elemental bismuth [5] 
to the first report of the quantum Hall effect at room temperature in graphene [6j. In 
this Letter we extend the quantum transport study of such Dirac structures to the critical 
point of a change in the bulk Fermi surface topology. From a semiclassical point of view, 
it is understood that magnetoresistance (MR) is sensitive to the topology of given Fermi 
surfaces [Tj: the magnetic field B drives the electrons in orbits around the Fermi surface 
(FS), sensing its geometry and topology. We here describe the effect of the Dirac structure 
in this context of Fermi surface topology, which we observe as a magnetoresistance effect of 
pure quantum origin across the bulk Dirac node. 

We study the system BiTel that possesses a Dirac node generated by Rashba-type spin- 
orbit coupling ol ■ (<r x k) 0i, where a is the structure-specific Rashba parameter, <t and 
k are the spin and momentum operator, respectively. This layered semiconductor breaks 
inversion symmetry by crystallizing in the polar space group P3ml so that the spin-orbit 
interaction takes the form of a Rashba term in the bulk 3D band structure. This has 
been confirmed by angle-resolved photoemission spectroscopy |10j and relativistic ab inito 
calculations nn- In this system, both the conduction band minimum and the valence band 
maximum are located near the A point in the hexagonal prism-shaped Brillioun zone (Fig. 
1(b)); taking A as the origin, the conduction electrons near the band edge can be described 
by the Hamiltonian 
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where the momentum k is decomposed into kit in the A-L-H plane and k z parallel to A — T as 
well as a. The quasi-two-dimensionality of the crystal structure leaves "H dominated by the 
in-plane components. Typical band parameters for BiTel are mo = 0.09m e , |c»£| = 3.85eV • A 
and m z /m 0 ~ 5, as reported in photoemission [TO], optical [12] and transport [T3] studies. 
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FIG. 1. (Color Online) (a) Schematic evolution of BiTel Fermi surface (FS) from a spindle-torus 
(Ep > 0) through horn-torus (Ep = 0) to a ring-torus (Ep < 0). Dashed lines define the IFS. (b) 
Location of FS in the Brillouin zone, (c) Representative p(T ) of samples with Ep above (green), 
near (black) and below (blue) the Dirac point with Hall carrier density denoted in parenthesis. 

To emphasize the competition between the spin-splitting a. ■ (<t x fc||) and kinetic energy 
h 2 k^/2mo, Ep is zeroed at the band crossing point at A, which is also the neutrality point of 
an effective Dirac Fermion with the Fermi velocity Vp m 5.35 x 10 5 m/s m Here the Dirac 
point defines a change of Fermi surface topology (Fig. 1(a)) from a spindle-torus (Ep > 0) 
through a horn-torus (Ep = 0) to a ring-torus (Ep < 0). The inner Fermi surface (IFS) and 
outer Fermi surface (OFS) are defined at each k z slice as shown in Fig. 1(a); the OFS is 
always electron-like while the IFS is bipolar depending on Ep and k z . 

The simple band structure described by Eq. (JT|) and degenerate semiconductor nature 
(usually n-type) of BiTel allows a systematic exploration of this system upon doping, which 
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readily transforms the FS progressively as depicted in Fig. 1(a). Single crystals of BiTel 
were grown in a vertical Bridgman furnace, intentionally doped with Cu to improve electron 
mobility /i e and vary the carrier density n e (2]. MR and Hall effects are measured in four- 
probe configurations in a Quantum Design Physical Property Measurement System with B 
applied along the polar c-axis. Compared with preceding studies jT3l fT5] our BiTel samples 
exhibit relatively large p 300K /p 21v = 4 ~ 5.8 (Fig. 1(c)) and enhanced low temperature /j, e 
typically between 800 ~ 3000 cm 2 /V-s. n e varies within the range of 1.2 ~ 6 x 10 19 /cm 3 
for samples taken from the same ampoule, implying a composition gradient intrinsic to 
Bridgman growth [IT] . The microscopic role of Cu in improving the electronic quality 
of crystals is still being investigated. At the lowest temperature, we observe clear SdH 
oscillations over an extended range of Ep as shown for selected samples in Fig. 2(a). 

SdH oscillations is a common means to probe FS in metals m whose frequency / gives 
the extremal FS size S'*. via / = USk/2'Ke. In Fig. 2(a) SdH of two different / are traceable 
in all normalized MR curves (A p/p 0 = ( p(B ) — p(0))/p(0)): one consists of broad bumps 
starting from R ~ IT, the other oscillates rapidly for B > 10 T. These represent two 
coexisting FS extremas of distinct size and are assigned to the IFS and OFS at k z = 0 in 
accordance with [I314I5] . Ep is estimated [18] by comparing f° hS and / lts to Eq. (|T|) with 
mo = 0.095 m e , |a| = 3.85 eV-A, which are determined below. Following the descending 
f°FS f rom 414 q 1 ( sam pi e A) to 233 T (sample H), i.e., E F from 79 meV to -46 meV (see 
Fig. 2(c)), the held corresponding to the last IFS oscillation (defined as IFS quantum limit 
R^ L , denoted by triangles in Fig. 2(a)) approaches B — 0 in samples A to E, then moves 
toward higher fields in samples F to H, consistent with the closing and opening of the inner 
FS pocket across Ep = 0 in Fig. 1(a). 

Beyond the IFS quantum limit, we observe distinctly different behavior for Ep > 0 and 
Ep < 0. It is most illuminating to contrast in Fig.2 (a) samples E (18 meV) and F (—24 
meV) which bracket the Dirac point. In sample E, at B > B the MR grows rapidly above 
2 T with a steeper slope than B < B C ^ L , while in sample F, the MR at B > B^ h is strongly 
suppressed relative to B < R^ L . This distinction is systematically visible in samples B-D, 
G-H with higher B ^ L . Apparently the MR beyond the QL provides a sensitive probe of 
the sign of Ep and corresponding Fermi surface topology in BiTel. This quantum transport 
behavior and its connection to FS topology are the main result of this Letter. 

We suggest that the physical origin of this bifurcation is the quantum transport of the 
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FIG. 2. (Color Online) (a) Normalized MR of samples A-H. Curves are scaled and offset to 
emphasize Landau level-related patterns. Triangles denote the IFS quantum limit and the dashed 
lines are extrapolation of the MR below the quantum limit, (b) Simulated IFS density of states 
£)its f or sam pi es A-H. Curves are normalized and offset for clarity, (c) Ep of samples A-H relative 
to Dirac point, (d) In-plane dispersion at k z = 0 when B = 0 (upper axis) and selected Landau 
levels as a function of B (lower axis). The colors blue and green represent respectively the IFS and 
OFS Landau levels, (e) Magnified R-evolution of D IFS near the Dirac point. 


relevant Landau levels (LL) for the density of states (DOS) at Ep [f?]. Figure 2(d) shows 
the in-plane electronic states without (upper axis) and with (lower axis) B. The quantized 
LL energies are given by mm 


heB 
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where n is the LL index, including n — 0 (Eo = heB/2mo). Selected LLs (n = —40 ~ 4) 
are shown in Fig. 2(d) and are divided into two groups colored blue and green, whose 
successive intersection with Ep gives rise to the IFS and OFS oscillations, respectively. We 
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note that because E n (B) for n < 0 are non-monotonic, in increasing B these LLs experience 
a transition from fFS-like to OFS-like. We illustrate this with the change in color of the 
negative LLs in Fig. 2(d). Blue (IFS) LLs resemble that of a Dirac Fermion with a slight 
R-dispersive n = 0 LL due to the zero-point energy of the parabolic band with the effective 
mass mo. Green (OFS) LLs can be seen to restrict the overall chemical potential within a 
weak variation. 

In contrast to this quantization in transverse fc||, the longitudinal k z is unaffected, recalling 
a residual ID contribution from every LL satisfying E n < Ep |2Uj . As the OFS LLs are 
insensitive to the sign of Ep, we focus on the IFS to understand the observed MR behavior. 
The IFS density of states D lhS is given by: 
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with E n explicitly expressed in Eq. ([2]) and Ep taken as R-independent for each sample. 
In Eq. (J 3 J) , Y y/2m z /2'K%y/Ep — E n counts the states at each cyclotron motion guiding 
center and eB/2nh counts the number of guiding centers. Note that we include the partial 
contribution of the lowest LL consistent with our definition of IFS as discussed above. The 
simulated R Its for samples A-G from Eq. (J3| are displayed in Fig. 2(b) with each E n 
broadened with a Lorentzian 40 K wide. 

Figure 2(b) reproduces the essential results in Fig. 2(a) including the bifurcation between 
samples with Ep > 0 and Ep < 0. This remarkable agreement suggests that the MR 
behavior reflects a quantum Landau level effect and is a measure of D us . I 11 Eq. (J3| the 
inverse square-root 1 f\jEp — E n is most sensitive to the uppermost LL below Ep. Thus 
R IFS at B > R QL captures the difference between E 0 (Ep > 0) and R_ 1 (Ep < 0) as a 
function of B. In the vicinity of the Dirac point, Fig. 2(e) shows a magnified view of _D Its 
with B gradually turned on. At Ep > 0, the contribution to D lhS from the weakly dispersive 
n = 0 LL grows steadily, reflecting primarily the increasing capacity per LL proportional to 
B. At Ep < 0, the suppression of D lhS is due predominantly to the rapid -y/B dispersion 
of n — — 1 LL, which in the present B range overwhelms the R-lincar LL degeneracy term. 

The finite m z in the bulk Rashba band structure is key to this effect by acting as a 
particle-hole symmetry (PHS) breaking term. In the context of band structure, in addition 
to the deformation of the Dirac node due to m 0 , m z breaks PHS of the node by forcing both 
the in-plane electrons and holes to be electron-like in k z . This allows qualitatively different 
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MR properties of the Dirac node in BiTel at Ep > 0 and < 0, in contrast to the nodes 
in graphene [ 2 ] and Weyl semimetals [ 21 ], where Ep and —Ep yields identical MR profile 
due to PHS. In the context of quantized Landau levels, m z breaks PHS when Ep is exactly 
between two adjacent LLs of energies E a < Ep the inverse square root preferentially samples 
E a and thus enables measuring the dispersion of E a with B in transport, as is the present 
case. This manner of PHS breaking is general for Landau quantization in 3D systems m 
and has only recently been addressed in the context of the Nernst effect near the quantum 
limit in graphite [ 22 ] • 

Beyond dimensionality and symmetry, the interconnection by the spin with a large OFS 
also distinguishes the transport of the Dirac node in BiTel from that in graphene. Con¬ 
ventionally, it may be expected that the IFS and OFS would simply be additive in their 

•1 • 1 , . . /«-»■ -f^OFS -<-»-IFSk 1 A 1 T T r 

contributions to the conductivity tensor (a = a + a ) so that A a due to LL formation 
would follow D Its . The agreement of Fig. 2(a) and (b), however, shows empirically that 
Ap r\j D lhS . We hypothesize that this is caused by the large difference in size of the spin- 
polarized IFS and OFS, which becomes extreme near the Dirac point. Considering such a 
scattering phase space, for the IFS the inter IFS-OFS backscattering events dominate over 
the intraband backscattering. The major role of the IFS, among the total electrical current 
carried almost completely by the OFS, is then through D lhS modulations that affect the 
interband scattering rate 1 / ti_o (and thus p ). 

The above identification of FS topology is independently supported by the observed SdH 
oscillations, most unambiguously shown by comparing the FS radii k OFS and fc IFS at k z — 0 
as in Fig. 3(a). This captures the vanishing IFS at the Dirac point where / OFS — 320 T, as 
expected from the relationship between the two coupled Fermi surfaces: 

k OFS - k 1FS = 2k a for E F > 0 

(4) 

k OFS + k 1FS = 2 k a for Ep < 0. 

The term k a = mo\ct\/Ti 2 is the offset of band minimum generated by ol ■ (cr x fc||), which 
performs a pure translation (gray arrow in Fig. 3(a) inset) between two spin branches on 
any fc-space cut through the k z axis. Eq. Q is a natural consequence of Eq. ([Tj), confirming 
that this expression effectively describes the electrons near Ep = 0. k a is estimated to be 
0.049 ± 0.008/A. This is consistent with m 0 = 0.095m e , |a| = 3.85 eV-A which we have 
used in the calculation above and is in reasonable agreement (within 6 %) with previous 
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FIG. 3. (Color Online) (a) OFS and IFS radii at k z = 0 (fcoFS and &ifs) for all measured samples 
fitted with Eq. The Rashba bands are shown in the inset, (b) IFS index plot of samples 

A,B,C,D,G and H. The inset shows the oscillatory resistivity after background subtraction. 

reports Enina, reassuring that the Cu-doping which improves the carrier mobility preserves 
the Rashba band structure. The IFS SdH senses both positive and negative low-index LLs 
across the Fermi surface topology change, as shown in the index plot Fig. 3(b). 

We next discuss the temperature (T) dependence of the observed quantum transport 
behaviors. The effect of T on SdH oscillations is a useful probe of the Landau levels and 
associated carrier effective mass and scattering EH- We have found that Kohler’s rule p], 
which states that A p/p 0 depends only on the Hall angle uj c t and thus B / p 0 , provides an 
incisive tool for isolating the LL contribution to transport across the QL. The Kohler’s plot 
of MR curves taken at various T of four different Fermi surfaces (samples A,C,F,G) are 
shown in Fig. 4. First, for B < B® h , the MR curves of each FS fall onto a single trace 
despite oscillatory deviation, implying that the scattering is hardly affected by 5 or T. 
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FIG. 4. (Color Online) (a)-(d) Temperature evolution of MR captured by the Kohler’s plot of 
samples A,C,F,G. The inset of (d) shows the positions of Ep at 14 T for each sample. 

Taking this collapsed curve as the background p Back onto which the SdH is superimposed 
assists the IFS index analysis in Fig. 3(b). We note that a previous work |1TJ adopts a 
similar procedure in studying SdH in BiTel. 

The validity of Kohler’s rule within the QL restricts the most plausible cause of the 
violation (at high fields beyond the QL in Fig. 4(b)-(d)) to quantum effects. We suggest 
that these violations are consistent with the unconventional Dirac LLs, in particular the 
n — 0 and n = —1. MR of samples C and F are notably enhanced with elevated T, 
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reflecting thermal excitation from the n — 0 LL both upward and downward (see inset of 
Fig. 4(d)). In sample G, although the IFS SdH appears weaker compared to C or A, the 
suppressed Kohler’s curve upon warming implies that the broad shoulder like feature in the 
MR is an intrinsic peak associated with the meeting of the coherent n — — 1 LL with the 
Ep. Additionally, the contrary T-evolution of samples C and G despite their similar B^ L 
demonstrates the asymmetry of quantum transport above and below the Dirac node. In 
this way it can be seen that Kohler’s scaling brings out the LL effects beyond the QL which 
would otherwise be overlooked. 

The uniqueness of BiTel lies in the particularly large value of a. compared to other 
Rashba systems maim. This strengthens the Dirac Fermion behaviors observed above and 
determines the large inter-LL scale: at 14 T E\ — Eq ~ 88 meV while Eq — E_\ ~ 71 meV. 
Furthermore, the Dirac point is located ~94 meV above the conduction band minimum, 
making Ep < 0 accessible at an appreciable doping level without entering the localization 
regime [23j. We note that this Ep < 0 regime in Rashba systems is of intense theoretical 
interest from the standpoint of enhanced superconducting instabilities [24], spin torque 
efficiency [26] and superconductivity with peculiar symmetry ra- Our present transport 
identification of the Ep < 0 FS topology thus paves the way for realizing the material host 
for these exciting proposals. 

In conclusion, we have established the quantum magnetotransport properties of the Dirac 
node in BiTel, which is markedly asymmetric for Ep > 0 and < 0. This effect is the result of 
a change in the bulk FS topology. Kohler’s scaling is shown to be an effective tool to analyze 
the unconventional low-index Landau levels. We predict the latter effect will be useful in 
the general study of band crossings in multi-band systems. 
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